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Abstract 

We study the evolution of the particle density, dn/drj at fixed 
T], with the beam rapidity Y in the framework of string percolation 
model. Our main results are: (i) The width of the "plateau" increases 
proportionally to Y, (ii) limiting fragmentation is violated, and (iii) 
the particle density, reduces to a step function. 

We work in the framework of the string percolation and one of the main 
results is the presence of an extended "plateau" with a length Ay propor- 
tional to Y{Ay ~ lAY), where Y is the beam rapidity in the center of mass 
system: Y = ln{-\f{s) / rup) . This result favors saturation models [1-4] with 
formation of longitudinal fields (flux tubes, effective strings) and naturally 
explains the presence of long range (pseudo-rapidity) correlations and the 
ridge phenomenon. See [5] for a discussion. Consequences of dominance 
of the "plateau" are the violation of limiting fragmentation and evolution 
towards a step function. 

In the string percolation model multi-particle production in a two nuclei 
collision at high energy can be explained in terms of the produced strings 
along the collision axis, between the projectile and target. These strings de- 
cay into new ones by qq or qq- qq pair production and subsequently hadronize 
to produce the observed hadrons. Due to confinement, the color of these 
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strings is confined to a small area in transverse space Si — ttTq , with 
Tq ^ 0.2- 0.3 fm. With increasing energy and/or atomic number of the 
colliding particles, the number of strings Ng grows and they start to overlap 
forming clusters, very much like disks in two-dimensional percolation theory. 
At a certain critical density, a macroscopic cluster appears, which marks the 
percolation phase transition [6] [7] . 

In string percolation the relevant quantity is the transverse impact pa- 
rameter density r;* which, in the case of pp collisions, we write as 

- {^m, (1) 

where tq is the single string transverse size, 7?^ ~ 1 fm is the proton transverse 
size and iV* is the average number of single strings. For values of ?7* below the 
critical 2- dimensional density for percolation, 77* ~ 1.15 — 1.5 [8], the formed 
strings do not interact and collective effects are not present. For values 
1]^ > ?7* one observes the formation of long strings due to fusion, stretching 
between the beam and the target. We assume that 1]^ > if^. 

The particle density dn/dy at mid-rapidity is related to the average num- 
ber Ng of strings 

^ - F{v')Ns, (2) 

where ^(77*) is the color reduction factor [9], due to color summation of 
random colors, 

/ 1 _ p-v* 

For small 77*, F(r]*) 1, and for 77* > rji, Firf) —k=. Starting with an 
exponential growth of the average number Ng of strings, 

Ns ~ e^^^, (4) 
with A ~ 0.2 — 0.3 [10,11], we write for the particle density, at y ~ 0, 

and for the full rapidity distribution, [12,13], 

dn, wdn,. 
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where is the single string density, ^\pp is the pp density, and a is a 
constant, depending on the nature of produced particles. 

Note that the exponential behavior in (4), iV* ~ e^^-^, is not an assump- 
tion but results from a simple application of conservation of energy to dn/dy 
- making use of (2), (3) and the high density hmit - to obtain (4) with 
A ~ 2/7, [11]. The rise of the multiphcity plateau is not proportional to Y 
but to exp(AF). 

Following [12,13] we write, for rj >0, 

dn.,^ 1 

dy\ e^^^^ + l' 

a and 5 being free parameters. Finally, using the Jacobean J of the y ^ rj 
transformation, we construct dn/dr)\pp, 

dn dn , . 

, 2, 2 

with J = cos T] ^ ^ _ m ^ assumption fix k at the effective 

k+sinh^r] Pt 

value 1.2. It corresponds to m-jr ~ 0.14 GeV and ~ 0.3 GeV also in 
agreement with [14]. 

In Fig. 1 we show our fits to dnch/dt] for pp collisions (excluding single 
diffraction) from 53 GeV to 1.8 TeV [15], and at LHC energies [16]. The free 
parameters where fixed at A = 0.23 ±0.005, 6 = 0.61 ±0.15, a = 0.27 ±0.03, 
a = 0.8 ±0.2. 

In Fig. 2 we show the evolution with Y at different values of pseudo- 
rapidity, rj. Evolution within the "plateau" is slower than evolution in the 
"fragmentation" region. 

Making use of (6) and (7) it is easily seen what happens. If 77 > is in 
the "plateau" region dn/drj evolves as e^^ , 

(IT) 

^(77<(1-Q;)y,y)~e^^. (9) 

(JjI j 

If \'q\ is in the "fragmentation" region dn/dr] evolves as e^^"*"^^^ 

^{V>{1- a)Y,Y) ^ e^^+^^. (10) 
dr] 

For fixed r] with increasing Y, at some stage one falls in the situation (9), as 
seen in Fig. 2 . 
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Figure 1: dNch/drj at different ^/s. For data points see references [15], [16]. 
Data from UA5, 546 GeV and UAl at 540 GeV being inconsistent, were not 
included. Data from P238, 630 GeV are also not included. Note that the 
^/s = 2.36 TeV and 0.9 TeV plotted is from CMS for clarity of the plot the 
ALICE results also follow in these range and they differ form the CMS data 
for less than 0.2. 
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Figure 2: Comparison of the results from the evolution of the drich/dr] with 
the beam rapidity from equation (7) for different fixed rj values. Blue circles 
are used for data at r/ = 2, red squares are used for data a.t rj = 2.5, and 
green triangles are used for data at 1] = 3 [15, 16]. 
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Our main result is that if we start at a given Yq with ^|r?>o < ^\v=o 
after enough Y evolution we end up with ^|r?>o — ^k=o (see Fig. 2, where 
we have worked not with ^ but with This means that a "plateau" 

proportional to Y is developing. 

Such "plateau" may be responsible for forward- backward long range 
rapidity correlations and the ridge phenomenon. It is well known, see for 
instance [17,5], that in simple models the forward-backward correlation pa- 
rameter b depends on the "plateau" density. If we have a long distance true 
plateau, then b is constant over a long distance. On the other hand the height 
of the ridge structure or the near ridge particles is also dependent on particle 
density. The same occurs in the frame of Color Glass Condensate model [18]. 

We discuss next the question of limiting fragmentation. For 77 large, 
77 > (1 — a)Y we have 

dn yY I (^-°) v V 
— ~ e s~'^ s 

dv (11) 

^ eA(y-„)+(i^r+(i^(-r,) 



The limiting fragmentation condition, i.e. the particle density being only a 
function of F — r;, is 1 — a = 1 — X6, or 

a/X = 5. (12) 

Experimentally a/X > I and S < 1 and including errors, a/X = 1.17±0.15 
and 6 — .61± 0.15, we see that (12) is not satisfied. 
To see what is expected we write (11) in the form 

^ ^ e^(y-r>)+^iy-n)+[^-^]{-v)^ (13) 
dr) 

and ^ — 7> as (or Y) goes to infinity. In this limit the distribution becomes 
close to a step function. It looks as if the fast particles in the fragmentation 
region disappear to feed up the front region in the plateau. 

The arguments developed here for pp scattering, also apply to AA and 
pA collisions, to the extent that they have similar longitudinal structure. 
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